Phase sensitive measurements of order parameters for ultracold atoms 
through two particles interferometry 
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Nontrivial symmetry of order parameters is crucial in some of the most interesting quantum many- 
body states of ultracold atoms and condensed matter systems. Examples in cold atoms include 
p-wave Feshbach molecules and d-wave paired states of fermions that could be realized in optical 
lattices in the Hubbard regime. Identifying these states in experiments requires measurements of the 
relative phase of different components of the entangled pair wavefunction. We propose and discuss 
two schemes for such phase sensitive measurements, based on two-particle interference revealed 
in atom-atom or atomic density correlations. Our schemes can also be used for relative phase 
measurements for non-trivial particle-hole order parameters, such as d-density wave order. 



The concept of order parameter, which charac- 
terizes states with spontaneously broken symmetries, 
has been successfully applied to a wide range of 
physical phenomena such as the Higgs mechanism 
in high energy physics [lj], superfluidity in neutron 
stars0|, superconductivity Q, gaseous Bose-Einstein 
condensates [J] and charge and spin ordering in electron 
systems @. Recent works on condensed matter systems 
emphasized that order parameters can often be char- 
acterized by non-trivial orbital symmetries. For exam- 
ple, in contrast to conventional superconductors, which 
have isotropic s-wave electron pairing, high Tc cuprates 
exhibit d-wave pairingsQ, while superfluidity of 3 He 
or superconductivity in Sr2Ru04 exhibit triplet p-wave 
pairings;?]. Other examples of order parameters with 
non-trivial orbital symmetries discussed in the literature 
are high angular momentum Pomeranchuk instabilities of 
electron systems [1] and unconventional charge and spin 
density wave states Despite of the interests of such ex- 
otic states, the experimental verifications of these states 
are yet a challenging problem. Standard thermodynamic 
and transport properties can be used to observe the mod- 
ulation in the magnitude of quasiparticle gaps,but not 
the change of the sign of the order parameters 6] . Only 
phase sensitive experiments, such as the observations of 
Josephson effects in corner SQUID junctionsflfj and n- 
ring tricrystal experiments [llj, have been considered as 
the definitive proof of the unconventional pairing for both 
cuprates and ruthenates 12]. In the case of states with 
anisotropic charge and spin orderings, the lack of exper- 
imental tests is one of the main reasons that their exis- 
tence remains an open question. 

During the last few years, a considerable progress 
has been achieved in creating analogues of strongly- 
correlated electron systems, using ultracold atoms in op- 
tical lattices (see refs. [l3| for reviews). One of the most 
challenging problems, which could be addressed in the fu- 
ture experiments, is the search for d-wave pairing in the 



repulsive Hubbard model[l4|. Realizations of other ex- 
otic states in cold-atom systems, such as <i-density wave 
states[15]. have been theoretically proposed. These states 
are characterized by order parameters with non-trivial 
angular dependence of the relative phase between the 
components of the entangled wavefunction. Hence, it 
is important to understand how tools of atomic physics 
can be used to perform tests of such quantum many-body 
states of ultracold atoms (l6|. 

In this paper, we discuss a scheme for performing such 
phase sensitive measurements. It is based on the analy- 
sis of atom-atom correlations resulting from two-particle 
interfere nce 11 7j . Our proposal builds on the theoreti- 
cal ideas [18| of using noise-correlations in atomic den- 
sity to characterize many-body states, and on the ex- 
perimental demonstration to measure atom-atom corre- 
lations, or atomic density noise spectroscopy with ultra- 
cold atoms[U|22|. This method should provide an un- 
ambiguous evidence for non-trivial pairings, including p- 
and (i-wave[l4. 23 . 24|, as well as for non-trivial particle- 
hole correlations such as in a d-density wave statefil Hij. 
It should also allow the direct observation of two particle 
coherence and nontrivial angular momentum of ultracold 
diatomic molecules. For example, for p-wave Feshbach 
molecules realized in JILA[25j], our approach should dis- 
tinguish between p x + ip y and p x — ip y states |24j|. 

We start by considering a Feshbach molecule, which 
consists of a pair of atoms, and has the center of mass 
momentum equal to zero 
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The two atoms making up the molecule can be either 
bosons or fermions. For concreteness, in this paper we 
focus on the case of two fermions in different hyperfine 
states labeled by a =tl, in analogy with states of a spin 
1/2 particle. Here tp(k) is the wavefunction of a molecule, 
c k<r i s a creation operator of a fermion atom in the state 
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FIG. 1: Illustrations of using two-atom interference to mea- 
sure the relative phase between different components of 
molecules after dissociation, a) Scheme I: Free propagating 
atoms are reflected in mirrors and mixed in beam splitters 
denoted by S. Coincidences are counted between detectors on 
opposite sides, e.g. Di and D3. b) Bragg pulses with wave- 
vectors G = p — q and — G are used to exchange (atomic mir- 
rors) or mix (atomic beam-splitters, S) components q "f and 
p t, as well as q 4. and p 4(|q| = |p|). In scheme II, the 
Bragg pulse applied at the beginning of expansion carries out 
reflections on mirrors and mixings in beam splitters in a single 
operation. 



with momentum k and hyperfme state a. The symme- 
try of ^(k) determines the nature of the paired state. 
We assume that the potential binding the two atoms 
is removed instantaneously and the released atoms sub- 
sequently evolve as free particles. Experimentally this 
can be achieved either by changing the magnetic field 
abruptly near a Feshbach resonance or by applying an 
RF pulse 2l|, |25j . The released pair of atoms are in a su- 
perposition of momentum (k, — k) pairs with amplitudes 
ip(k). Our goal is to find a method to measure the rela- 
tive phases between ^(k) for different k. 

Scheme I. We first explain the main idea through 
the scheme of Fig[T^, analogous to the quantum optics 
scheme of 26]. Atomic mirrors and beam splitters are 
used to reflect and mix states with momenta p and q on 
one side, -p and -q on the other side. Time and space re- 
solved detectors in opposite sides (e.g. Dl and D3) allow 
measurements of correlations resulting from the interfer- 
ence between ip(p) and ip(q), and thus, can reveal the 
relative phase between these components. The atomic 
mirrors and beam splitters are based on Bragg diffrac- 
tions from laser beams which make standing waves with 
wave vectors ±G, and they couple states with the same 
spin and magnitude of momenta (|q| = |p|) but whose 
momenta differ by p — q = ±G. For simplicity, we con- 
sider long Bragg pulses where a perfect Bragg diffrac- 
tion can be achieved, i.e. no other diffraction order is 
involved. Later, we demonstrate that the phase sensi- 
tive measurements can also be carried out with short and 
strong Bragg pulses which can introduce higher diffrac- 
tion orders. The atomic mirrors in Fig JT] can be realized 
through Bragg pulses whose amplitudes and durations 
are chosen to produce ir pulse so that it converts ±p into 
±q and vise versa. On the other hand, a n/2 pulse in- 



duces mixing between states with momenta ±p and ±q, 
realizing a beam splitter (denoted by S in Fig|T]). We 
can express the original fermion operators in terms of 
the operators after the mixing as follows: 

e~ iB ^cL. = cos/?/} - isin/3e ix t 4, 

e ~ iepT c pt = _i sin / 3e_iXT 4 + cos p 4 > 
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Here d\ are creation operators for particles observed in 
detectors (i = 1,...,4). The mixing angle /3 (of the 
order of tt/2) and spin-dependent phases \a can be con- 
trolled through the amplitudes, durations and relative 
phases of the Bragg laser pulses. We denote by 8^ the 
phase accumulated by an atomic component with mo- 
mentum k and spin a during the propagation between 
the source and the beam splitters. 

If we assume that molecular wavefunctions for wave 
vectors q and p differ only in phase, i.e. ip(k) = \ip\ e 1 ^, 
we find the following expressions for the coincidence 
counts of rii — d\di 
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and similarly for ( n% n 3 ) c and {n% 114, ) c . The oscillatory 
behavior of the correlation as a function of $/ probes the 
coherence of pairing in the molecule. To vary <£>/, one can, 
for instance, change the phases Xa- Moreover, if we know 
the precise value of such coincidence signals yield the 



relative phase 



between different molecular com- 



ponents. In the absence of precise knowledge of <!>/, the 
phase difference q — <p p could be extracted through a 
scheme analogous to white light fringes in classical op- 
tics, whose pattern and shape can reveal the existence 
of fundamental phase factors 27j. Note, however, that 
k dependence of the phase factors acquired during the 
propagation and the reflection may render these meth- 
ods unreliable. Thus, we consider a second scheme which 
avoids such a problem. 

Scheme II. In this alternative scheme, we apply a 7r/2 
Bragg pulse at the very beginning of the expansion to 
mix atomic components with momenta q t and p as 
well as — q 1 and — p 4- This realizes, in a single opera- 
tion, reflections on the mirrors and mixing on the beam 
splitters. In scheme II, there is a common mode prop- 
agation after the Bragg pulse, and phases acquired dur- 
ing the expansion do not affect interference. Two atom 
interference is revealed by coincidence counts with point 
detectors just as in the previous scheme. The scheme can 
be generalized to many-body case by replacing coincident 
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counts between point detectors with density imaging and 
studying noise correlations between patterns registered 
on opposite sides (see below) . 

To discuss scheme II, we start again with the exam- 
ple of a dissociated Feshbach molecule, described by the 
wavefunction in Eq.([T]). We consider the case in which 
the Bragg pulses for spin up and down atoms differ only 
in the phase, and such pulses are created by the potentials 
V CT (r) = 2Vb cos(G.r — Xa)- Here we assume again a per- 
fect Bragg diffraction. Detectors (i = 1, 2, 3, 4) detect 
atoms with momenta and spins p t, q ti ~°± b ~P 4-i re_ 
spectively. The only difference between this scheme and 
scheme I is the absence of phase factors e 10 ' 1 " . As a result, 
coincidence counts have forms similar to Eq.Q, with <f>/ 
replaced by <!>// = \t ~~ Xi- Therefore, provided that we 
know the phase difference Xi ~ Xi associated with the 
two Bragg pulses, atom-atom coincidence counts directly 
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the pairing symmetry in Eq.(fT|). 



Many-body state analysis. We now apply scheme II to 
a BCS state of fermions |\&) = Y\ k ( M k + w kcL.cl k ,) |0)- 
This BCS wavefunction is general and can describe 
weakly-coupled BCS paired states as well as a con- 
densate of tightly-bound molecules. Here, we con- 
sider the generic diffraction pulse that can mix states 
whose momenta are separated by any integer multiple 
of G. The effect of the mixing pulse is described by 
the transformation of particle creation operators: cj^ — > 

c kt — L~>m a 0,m t c k+mGt' C k-Gt ^ L k-Gt — 
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ak t,m e_l(m+1)xTc ik+mGt' and analogously for c [ _ ki 
and cL k+G ^. The scattering amplitudes o^" m are con- 
trolled by the diffraction pulse amplitude Vq and its du- 
ration t. We assume that before the mixing pulse, only 
states with momenta ±k, ±(k — G), which are close to 
the Fermi surface, have finite probabilities to be occu- 
pied, while states with momenta ±(k — mG) for m ^ 0, 1 
, which are far from the Fermi surface, are empty. The 
mixing pulse then induces the interference between par- 
ticles with momenta ±k and ±(k — G). 

The signature of non-trivial pairing of the BCS wave- 
function shows up in the angular dependence of the phase 
4>k in «k = l^kle 1 ^. In order to probe the relative phase 
A(f> = </>k — </>k-G between pairs with momenta k and 
k — G, we consider the following density noise correla- 
tion after the interference: 
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In analogy with the case of a Feshbach molecule, the first 
line in the RHS of Eq.(j4]) contains an interference term 
which depends on the relative phase A<f> as well as on 
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FIG. 2: a) In order to take into account the finite resolu- 
tion of detectors and the integration in absorption imaging, 
the density noise correlation is integrated over the cylinders 
shown in the figure, b) Integrated density noise correlations 
{SNv t SNvi) as a function of phase difference xt — Xi f° r a 
strong diffraction pulse of amplitude Vo/Er = 2 and a du- 
ration t which yields the maximum oscillation of the signal. 
Blue, Green(dash-dotted line) and Red(dashed line) curves 
correspond to A<f> = <f>^ — <^k-G = 0, 7r/2 and 7T, respectively. 



Space- and time-resolved single atom detection [li], 22, 
[28[ permits direct measurements of atom-atom correla- 
tions for specific momenta, corresponding to Eq.Q. Al- 
ternatively, one may look for noise correlation in absorp- 
tion images after time of flight [l8j. In this case, absorp- 
tion imaging, as well as finite resolution of detectors, re- 
sult in the integration of the atomic density. In order to 
take into account these effects, we have integrated Eq. (f5|) 
over ranges of momenta as shown in Fig|2^,. We present 
in Fig[2p the numerical result of this integration, which 
displays noise correlation in integrated density vs. the 
phase difference Xt — Xi °f the diffraction pulses. Here 
we took the integration range to be \Ak y \ = |G|/10, 
|Afc x | = |G|/10, \Ak z \ = 5|G| and the pairing gap to 
be A w O.lEp. The diffraction pulse amplitude is set to 
Vo/Er = 2 where En — |G| 2 /8m is the recoil energy, 
and its duration is chosen to have the maximum oscilla- 
tion of the signal. We assume that the integration range 
is sufficiently small that the phases of the Cooper pairs 
<pk and 0k-G are constant in the integration range. 

The oscillatory behavior of the integrated noise cor- 
relations (SNv r 6Nv±) as a function of xt ~ X4.(See Fig 
[5p) should provide an unambiguous proof of the Cooper 
pair coherence. Moreover, the value of the correlation 
at Xi ~ Xi = yields information about the phase dif- 
ference Acf> — (j)^ — </>k_G, which is the quantity we are 
interested in. The value of the correlation at Xi ~ Xi = 
also depends on the scattering amplitudes oq a m , and thus 
on Vqt. In Fig[3l we present the integrated noise corre- 
lation signal (5Ny r SNy i } at xt — Xi = as a function 
of Vqt for three different values of Acf>, and find striking 
differences. We conclude that it should be possible to 
discriminate between A(f> = and A<j) = it even when 
full 3D resolution is not available. 

In discussions so far, we assumed that the BCS pairs or 
molecules are at rest before dissociation. When molecules 
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FIG. 3: Integrated noise correlations (<5AV t 5Ny±) as a func- 
tion of VoT a) for a Bragg pulse amplitude Vq/Er — 2, 
and b) for a Bragg pulse amplitude Vo/Er = 20. Blue, 
Green (dash-dotted line) and Red(dashed line) curves corre- 
spond to A(f> = 4>k — 0k-G = 0, 7r/2 and -k, respectively. 




FIG. 4: a) Illustration of the sign of the CDW amplitude 
^ph(k) = (cj^Ck+Qo-) for a!- wave CDW state, b), c) Phase 
sensitive detection of the symmetry of the rf-wave CDW state 
in TOF experiments. For b), two pulses which transfer mo- 
menta G and G' are applied at the beginning of expansion. 
In c), a single pulse with momentum transfer G is applied. 
All the couplings through the Bragg pulses are indicated by 
solid arrows. Here Q is the wave vector of CDW. 



are cold but not condensed, there is a spread in the center 
of mass momenta determined by the temperature. Even 
in this case, there is still a perfect coherence between dif- 
ferent parts of the wavefunction of each molecule, yield- 
ing a two-body interference. However, the average of 
these interference terms over the center of mass momenta 
of individual molecules could potentially result in the 
washing out of noise correlations. We expect this sup- 
pression to be moderate as long as the average thermal 
momenta of molecules is smaller than the characteristic 
momenta of expanding atoms. 

Systems with particle-hole correlations. There are sev- 
eral types of many-body states characterized by corre- 
lations in the particle-hole channel, such as charge and 
spin density wave states (CDW and SDW). The most 
exotic of them have a finite angular momentum. This 
means that we have (c ko .Ck+Qo-) = if> p h(k), where -0pft,(k) 
has a non-trivial angular dependence. Our scheme above 
can be generalized to provide an unambiguous phase 
sensitive detection of such states as well. To be con- 
crete, let us consider a 2D system near half-filling. In 
this case, one can combine two different measurements 
of correlation functions to obtain the information on 
the order parameter ip p h(k), as shown in FigH]D and 
c. In FigJ3t>, two Bragg pulses couple k and k' + Q, 



as well as k' and k + Q. Here, the correlation function 
(SrikSrik') contains an interference term proportional to 
ip p h(k)ip p h(]c.'). In FigHJ:, a Bragg pulse couples k and 
k', and the correlation function (SrikSn^i + q) contains the 
term i/)* ;i (k)^p^(k'). When combined, these information 
should not only provide evidence of the angular depen- 
dence of CDW, but also allow one to distinguish site and 
band centered density wave states. 

Conclusion. In this paper, we have proposed a new 
method for performing phase sensitive measurements 
of non-trivial order parameters in systems of ultra-cold 
atoms, with a view toward studying open problems in 
strongly correlated systems. Note that in contrast to 
scheme I, which was introduced in analogy to a scheme 
first introduced in Photon Quantum Optics, scheme II is 
specific to Quantum Atom Optics, and takes advantage 
of the unique features of cold atom systems. 

We acknowledge useful discussions with E. Altman, I. 
Bloch, M. Lukin, and thank C.I. Westbrook for his sug- 
gestions on the manuscript. This work was supported 
by the NSF grant DMR-0705472, Harvard MIT CUA, 
DARPA OLE program, AFOSR MURI, CNRS, ANR, 
Triangle de la Physique, IFRAF. 



[I] S. Weinberg, The Quantum Theory of Fields , (Cam- 
bridge University Press, 1996). 

[2] G. Baym and C. Pethick, Annu. Rev. Nucl. Sci. 25, 27 
(1975). 

[3] R. Schrieffer Theory of Superconductivity, (Perseus 

Books, Massachusetts, 1983). 
[4] F. Dalfovo, S. Giorgini, L. Pitaevskii, and S. Stringari, 

Rev. Mod. Phys. 71, 463 (1999). 
[5] O. Zachar, S. A. Kivelson, and V. J. Emery, Phys. Rev. 

B 57, 1422 (1998). 
[6] D. J. Scalapino, Phys. Rep. 250, 329 (1995); C. C. Tsuei 

and J. R. Kirtley, Rev. Mod. Phys. 72, 969 (2000). 
[7] A. Leggett, Rev. Mod. Phys. 47, 331 (1975); A. Macken- 
zie and Y. Maeno, Rev. Mod. Phys. 75, 657 (2003). 
[8] V. Oganesyan, S. A. Kivelson, and E. Fradkin, Phys. 

Rev. B 64 195109 (2001); C. J. Wu and S. C. Zhang, 

Phys. Rev. Lett. 93, 036403 (2004); C. J. Wu, K. Sun, E. 

Fradkin and S.C. Zhang, Phys. Rev. B 75 115103 (2007). 
[9] S. Chakravarty et al, Phys. Rev. B 63, 094503 (2001); 

P.A. Lee, N. Nagaosa and X. G. Wen, Rev. Mod. Phys. 

78, 17 (2006). 

[10] D. A. Wollman et al, Phys. Rev. Lett. 74, 797 (1995). 
[11] C. C. Tsuei et al, Phys. Rev. Lett. 73, 593 (1994). 
[12] K. D. Nelson et al, Science 306, 1151 (2004). 
[13] I. Bloch, J. Dalibard and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008); M. Lewenstein et al, Advances in Physics 

56 Nos. 1-2, 243 (2007). 
[14] W. Hofstetter et al, Phys. Rev. Lett. 89, 220407 (2002). 
[15] C. Honerkamp and W. Hofstetter, Phys. Rev. Lett. 92, 

170403 (2004). 

[16] V. Gritsev, E. Demler and A. Polkovnikov, Phys. Rev. 
A 78, 063624 (2008); I. Carusotto and Y. Castin, Phys. 
Rev. Lett. 94, 223202 (2005). 



■5 



[17] M. O. Scully and M.S.Zubairy, Quantum Optics, (Cam- 
bridge University Press, 1997); N.R.Thomas et al, Phys. 
Rev. Lett. 93, 173201 (2004); Ch. Buggle, et al.,¥hys. 
Rev. Lett. 93, 173202 (2004). 

[18] E. Altman, E. Demler and M. D.Lukin, Phys. Rev. A 70, 
013603 (2004). 

[19] M. Schellekens et al, Science 10, 648 (2005); T. Jeltes et 

al, Nature 445, 402 (2007); M. Yasuda, and F. Shimizu, 

Phys. Rev. Lett. 77, 3090 (1996). 
[20] S. Foelling et al, Nature 434, 481-484 (2005); T. Rom 

et al, Nature, 444, 733-736 (2006); I. B. Spielman, W. 

D. Phillips and J. V. Porto, Phys. Rev. Lett. 98, 080404 

(2007) ; V. Guarrera et al, Phys. Rev. Lett. 100, 250403 

(2008) . 

[21] M. Greiner et al. Phys. Rev. Lett. 94, 110401 (2005). 



[22] A. Perrin et al, Phys. Rev. Lett. 99, 150405 (2007); K. 

Molmer et al, Phys. Rev. A 77, 033601 (2008). 
[23] D. W. Wang, M. D.Lukin and E. Demler, Phys. Rev. A 

72, 051604(R) (2005). 
[24] V. Gurarie and L. Radzihovsky, Annals of Physics 322, 

2 (2007). 

[25] J.P. Gaebler et al, Phys. Rev. Lett. 98, 200403 (2007). 
[26] M. A. Home, A. Shimony and A. Zollinger, Phys. Rev. 

Lett. 62, 2209 (1989); J. G. Rarity and P. R. Tapster, 

Phys. Rev. Lett. 64, 2495 (1990). 
[27] M. Born and E. Wolf, Principles of Optics, (Cambridge 

University Press, 1999). 
[28] R. Bucker et al, New Journal of Physics 11, 25 (2009). 



